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Probing magnetic order in ultracold lattice gases 
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A forthcoming challenge in ultracold lattice gases is the simulation of quantum magnetism. That 
involves both, the preparation of the lattice atomic gas in the desired spin state as well as the 
probing of the state. Here we demonstrate how a probing scheme, based on atom-light interfaces, 
gives access to the order parameters of non trivial quantum magnetic phases, allowing to characterize 
univocally strongly correlated magnetic systems produced in ultracold gases. The method, which 
furthermore is non demolishing, yields spatially resolved spin correlations and can be applied to 
bosons or fermions. As a proof of principle, we apply the method to detect the complete phase 
diagram displayed by a chain of (rotationally invariant) spin-1 bosons. 

PACS numbers: 37.10. Jk, 75.10.Pq, 42.50.-p 
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The possibility of using ultracold gases as quantum 
simulators, i.e. an experimental system that mimics 
models of condensed matter, high energy physics or quan- 
tum chemistry, was boosted by the seminal realization 
of Bose-Einstein condensate in alkaline atoms [1[ . Since 
then, the field has progressed enormously, both theo- 
retically and experimentally, moving from the weakly- 
correlated regime towards the strongly-correlated one 
[3, y|- This has been accomplished mainly by the use 
of optical lattices and Feshbach resonances. A forthcom- 
ing step is the simulation of quantum magnetism -spin 
ordering- in ultracold lattice gases. Recent experiments, 
showing super-exchange interactions in two nearby sites 
or novel methods to lower the temperatures enough to 
achieve this goal have already been reported [4j. 

One of the major difficulties in simulating magnetic 
systems is not only the preparation but also the faithful 
detection of the many-body state. A tomographic re- 
construction of the state is normally unrealistic due to 
the number of particles involved in the system. However, 
in quantum magnetism, ordering is usually reflected by 
a set of order parameters. An order parameter is a lo- 
cal operator A with a non vanishing expectation value, 
(</?|^4|(/j) =^ 0, that accounts for some broken symmetry 
embedded in the ground state \<p). Landau introduced 
the order parameter concept as a way to quantify the 
dramatic transformation of matter at a (classical) phase 
transition and since then it has become a cornerstone 
concept in physics. Notice that there exist also quantum 
phases like e.g. topological, whose symmetries are not 
local but global and, therefore, cannot be described by 
local order parameters. In these systems, strong quan- 
tum fluctuations prevents the formation of local ordering. 
While the total magnetization will be the appropriate 
order parameter for a ferromagnetic state ((M z ) =/= 0; 
M z = J2 i S z i, where S z i is the operator of the ith spin), 
the characterization of spin ordering in systems with zero 
magnetization requires different order parameters. Some- 
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FIG. 1: (Color online) Schematic detection setup, atoms 
placed in an optical lattice of periodicity d/2 (not shown) are 
illuminated by a laser beam in a standing wave configuration 
(yellow region) shifted by a from the optical lattice configu- 
ration. The output light is redirected through a polarimeter 
which measures its polarization through a homodyne detec- 
tion (HD). 



times order parameters are not observables (e.g. the su- 
pcrfluid order parameter) and cannot be directly mea- 
sured. Notice that any local (or global) order parameter 
could be extracted from the knowledge of all order spin- 
spin correlations since the later is equivalent to knowing 
the ground state. Detection methods are, however, re- 
stricted to lower moments. For instance, standard detec- 
tion methods like noise-correlations [5|, 1|| detect density- 
density correlations which can be linked to spin-spin cor- 
relations, but are not suitable for extracting higher order 
spin correlations. As proposed in pj, Bragg scattering 
of light gives rise to the spin structure factor, useful to 
detect Neel antiferromagnetic states, but will fail to de- 
tect other states that have zero magnetization and do not 
break time-reversal symmetry. An alternative technique 
for imaging spin states in optical lattices has been put 
forward [8(. Single lattice detectors |9| should be able 
to resolve magnetic ordering although they are presently 



limited to the type of atoms they can detect and consti- 
tute a destructive technique. 

A novel tool to detect magnetic correlations in lat- 
tice gases, based on quantum polarization spectroscopy 
(QPS), can overcome these limitations. The method, re- 
ported in Ref. [l0( and depicted schematically in Fig. [TJ 
takes advantage of the Faraday effect occurring when 
off-resonant light interacts with an atomic system (ei- 
ther bosonic or fermionic) with internal spin degrees of 
freedom. A polarized light pulse sent to a spin-polarized 
atomic sample couples its polarization to the atomic spin. 
As a result of the interaction, the light polarization ro- 
tates proportionally to the sample magnetization. In- 
deed, QPS is based in an effective dipole interaction 
Hamiltonian where an adiabatic elimination of the ex- 



cited states has been performed [11 
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The coupling constant k depends on the optical depth 
of the atomic sample and on the spontaneous emission 
rate. A light pulse propagating along the z-axis and 
strongly polarized in the x-direction is conveniently rep- 
resented by Stokes operators si = ^(a x a x — ala y ),S2 = 
\{ala x + a\.ay),S3 = ^{a\a x -a x a y ) where a x (a y ) is the 
annihilation operator of a photon with x (y) polarization. 
The atomic magnetic moments are included in the effec- 
tive angular momentum: J° ff = l/vL^2 n c n S zn where, 
L is the total number of spins . For strongly localized 
atomic Wannier functions, the coefficients c„, related to 
the standing wave intensity profile, can be approximated 
by: c„ = 2cos 2 [kpd(n — a)} with kp being the wave 
number of the probing laser, d the lattice spacing, and 
a the dimensionlcss shift between the probing laser and 
the lattice. The operator J° is mapped onto the light 
quadrature X = J dts 2 j \/N p h where we have assumed 
the incoming pulse to be strongly polarized along the x- 
direction, (( f dtsi)) — N p h ^ 1. The result of the map- 
ping is JTl| : 
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where X out (Xi n ) is the quadrature of the outgoing (in- 
coming) light. A homodyne measure of X out gives access 
not only to the effective angular momentum J° but also 
to the variance (as well as higher moments), demonstrat- 
ing the inherent quantum character of the probe. Let us 
focus on the variance e(kp, a) to prove how QPS can be 
linked to different order parameters: 



cff\2 



s(k P ,a) = (AJf) 



L^ 



c m c n Q z (m,n). (3) 



where Q z {m,n) = (S zm S zn ) - (S zm )(S zn ) is the two- 
point spin correlation function. As noticed in 12j, 
the magnetic structure factor, defined as S(q) = 
l/L J2inn cxp[iqd(m — n)](S zm S zn ), can be connected to 



the variance s(kp, a) as the average over a: 

e(k P ) = J dae(k P ,a) = \s(2k P ). (4) 

Although the magnetic structure factor contains the two- 
body correlation functions, it does not provides the same 
information as order parameters since by averaging over 
the phase shift a one loses some information on the cor- 
relations. In this work instead we assume an accurate 
control on the shift a and we introduce the quantity: 



Ae(k P , ai,a 2 ) = e(k P ,ai) - e(k P , a 2 ) 



(5) 



which is the difference of the signal with fixed wavevector 
kp at two different phase shifts. By appropriately choos- 
ing the parameters kp , u\ , a 2 , this quantity is linked to 
local order parameters able to detect not only ferro and 
antiferromagnetic (Neel) order, but also more exotic spin 
ordering. This makes QPS a valuable tool for quantum 
simulators. We support our claim analyzing a one dimen- 
sional spin-1 chain, a model characterized by different 
magnetic orderings and that can be realized with ultra- 
cold atoms in optical lattices. 

Spin-1 atoms confined in a deep optical lattice are well 
described by the Bose-Hubbard Hamiltonian [13}. For 
unit filling and for sufficiently small tunneling the sys- 
tem is in a Mott insulator state with one atom per site. 
Virtual tunneling of the atoms between neighboring sites 
gives rise to an effective magnetic interaction described 
by the bilinear-biquadratic Hamiltonian 1 1 31 ] : 



Hbb =J^cos(e)Si ■ S i+1 + sm(e)(Si ■ S i+1 ] 



(G) 



Hamiltonian ([6]) is derived within second order perturba- 
tion theory in the tunneling rate and the angle 9 e [— 7r; 7r] 
depends only on the two-body short range atom-atom in- 
teraction. This model has a rich phase diagram depend- 
ing on the angle 9 and has been extensively studied in the 



literature, see 14-181 and references therein. We review 



the phase diagram, depicted in Fig. HJa): 

The ferromagnetic phase.- For tt/2 < 9 < 57r/4, the 
ground state is ferromagnetic exhibiting a spontaneous 
magnetization, which serves as a local order parameter. 
For the remaining values of 9, the ground state has no 
net magnetization but different quantum phases: 

The critical phase.- In the interval ir/4 < 9 < tt/2 
the system is gapless due to soft collective modes at mo- 
menta q = 0,±27r/(3(i). The ground state spin-spin cor- 
relation functions (S Z iS z ^-\- r )) show period-3 oscillations 
[18j . In momentum space, this feature emerges as a peak 
at q = 2ir/(3d) in the mag netic structure factor S(q). Re- 
cently Lauchli et al. [15( have shown that nematic (i.e. 
quadrupolar) correlations at momentum q = 2ir/(3d) are 
enhanced in the critical phase while spin correlations be- 
come smaller when increasing 9 from 0.2ir to 0.57T. To- 
gether with the absence of the gap, the enhanced nematic 
correlations are a distinctive feature of the critical phase. 




FIG. 2: (Color online) (a) Phase diagram of the model of 
Eq.©. (b) The quantity e(7r/3d, a) (in units of k 2 ) for 9 — 
0.3-7T for L — 132 as a function of a. 



The Haldanc phase.- The interval — 7r/4 < 9 < 7r/4 
contains for 9 = the spin-1 Heisenberg chain and for 
tan(0) = 1/3 the exactly solvable AKLT point 0. The 
presence of a gap excludes local quasi long range order 
and spin correlations decay exponentially. The Haldane 
phase can be characterized by a hidden topological order 
parameter, called the string order parameter [19J, that 
cannot be revealed with local measurements. 

The dimcr phase.- In the region — 3ir/4 < 9 < — 7r/4 
the ground state is gapped and breaks translational in- 
variance, organizing in slightly correlated dimers. For 
— 37r/4 < 9 < — 7r/2 it is still under debate whether the 
system is always dimerized or it becomes nematic [20|. 
Numerical results [lj, [l5J show that the dimcr order pa- 
rameter, D = | (Hi — H i+ i) | where Hi = cos(9)Si ■ Si + i + 
sin(9)(Si ■ Si + i) 2 , is different from zero up to values very 
close to 9 = — 37r/4. 

Here we calculate the variance e(kp, a), which depends 
on all possible two-spin correlations for the ground state 
of the system at different points of the phase diagram 
given by 9. Ground states are found numerically using 
a density matrix renormalization group (see for exam- 
ple [2l|) with open boundary conditions. In Fig. [3] we 
show e(kp 1 a) for three different values of 9 correspond- 
ing to the critical, Haldane and dimer phases. A com- 
mon feature of the three phases is the presence of a high 
peak at kpd = ir/2 due to antifcrromagnetic correlations. 
Such peak is instead absent in the ferromagnetic phase. 
Apart from this, the plots in the three phases are quali- 
tatively different. In fact, in the critical phase, the signal 
is characterized by peaks at kpd ~ 7r/3 and kpd ~ 2ir/3. 
These resemble the peaks of the magnetic structure factor 
S(2kp) [ljj and are due to the period-3 oscillations of the 
correlation functions. These correlations are fundamen- 
tal for the detection of the critical phase. For 9 < — 7r/4 
we find the appearance of other small peaks at kpd = ir/4 
and kpd = 3-7r/4 signaling, as we show below, the dimer 
order. Finally, in the Haldane phase, the signal e(kp, a) 
lacks of non trivial features for kpd ^ n/2. 

We now use the quantity Ae(kp,a\, a^) introduced in 
Eq. (J5J) for the detection of the critical and dimer phase. 
To see how to choose the parameters kp,ai,(X2, let us 
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FIG. 3: (Color online) Left column, the function e(kp,a) (in 
units of ft 2 ) for different values of 9 in the three phases for 
L = 132: top 9 = -0.5tt (dimer), middle 9 = (Haldane), 
bottom 9 = 0.3-7T (critical). Right column, the same plots but 
restricted to a = 0. 



consider the critical phase. In this case it is natural 
to take kp = ir/Sd. Then we study the behavior of 
e(7t/3(i, a) in the critical phase as a function of a. 

The numerical analysis, reported in Fig. E^b), shows 
that the quantity e(7r/3d, a) is a sinusoidal oscillating 
function of a. The role of correlations at kp = n/3d is 
optimized by taking the difference between its maximum 
value at a± = 5/4 and its minimum at ai = 1/2. We 
then introduce C e = Ae(Tr/3d, 5/4, 1/2) which we con- 
sider as an order parameter for the critical phase. Since 
the ground state has zero angular momentum, we obtain: 
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The quantity C e is sensitive to correlations which oscil- 
late with a period 3 and represents a footprint of the 
critical phase. In Fig. |Ua) we show the signal C e for 
different values of 9. The results clearly show that the 
critical phase is very well detected by a positive value of 
C e . For 9 = 0.27T, in the Haldane phase and close to the 
phase transition, we still observe a large positive value, 
probably due to residual pcriod-3 correlations persisting 
in the Haldane phase. However, in this point, we find a 
non negligible dependence with the size of the sample. 
As shown in the inset of Fig. [UJa) a finite size scaling 
suggests that in the thermodynamical limit for L — > oo 
the quantity C e goes to zero as 1/L for 9 = 0.2ir, while 
for the other values of 9 > 0.247T it converges to a finite 
value. In Fig. H^a) we compare C e with the magnetic 




FIG. 4: (Color online) (a): The quantity C e = 
Ae(Tv/3d, 5/4, 1/2), in units of k 2 , (circles) as a function of 
6 for L — 132 and the rescaled magnetic structure factor 
[S(2n/3d) — l]/5 (squares). Inset: finite size scaling of C e for 
6 — 0.27T. The solid line scales as 1/L. (b): the quantity 
T> E = Ae(7r/4d, 1/2, 3/2), in units of k 2 , (circles) for different 
values of 6 for L = 132 and the rescaled dimer order parame- 
ter Dr = D/5.38. In the two plots, we distinguish the model 
phases with different shading: horizontal lines (dimer), no 
shading (Haldane), oblique lines (critical). 



structure factor S(2iv/3d). We confirm the results found 
in [15( | for chains up to L = 18: since the maximum is 
at 9 ~ 0.27T, the magnetic structure factor alone is not 
sufficient to distinguish the critical point at 6 = 0.257T. 
As shown in 15|, one would also need the quadrupolar 
structure factor to find the critical point. Our findings 



indicate that measuring only C e , which behaves as an or- 
der parameter, is sufficient to infer the occurrence of the 
phase transition. 

In the dimerized phase, the presence of peaks at kpd = 
it/A signals the pairing of neighboring spins. By averag- 
ing the signal e(kp,a) over a, these peaks disappear. 
Therefore, these features are not visible in the magnetic 
structure factor. Following a similar reasoning that leads 
to the definition of C e , we find that the quantity 
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is suitable for the detection of the dimer phase. The fac- 
tor sin[7r(m + n)/2] ensures that only the pairs of spins 
with positions m and n of opposite parity contribute to 
T> e . The differential signal V e is therefore an extension 
to long range correlations of the dimer order parameter 
D. In Fig. [4jb) we compare the signal T> e and the dimer 
order parameter D for different values of 6. Similar to 
D, the quantity T> 6 is significantly different from zero 
only in the dimerized phase, with the advantage of being 
experimentally measurable. 

The simulation of strongly correlated magnets with op- 
tical lattices can be foreseen in the next generation of 
experiments with ultracold atoms. Once this is accom- 
plished, the measurement of order parameters character- 
izing these phases will be of utmost importance. In this 
work we proposed the use of QPS for the direct observa- 
tion of magnetic order parameters in ultracold atoms in 
optical lattices. As a proof of principle, we have shown 



that this method allows us to unambiguously reconstruct 
the rotational invariant spin-1 chain phase diagram. We 
emphasize that the method is very general and could 
be applied to other spin models. Notice also that the 
QPS technique we propose is not based on any previous 
knowledge of the quantum phase to be characterized, but 
rather takes advantage of the features displayed by the 
probing scheme properly optimized. Although we are 
not able to measure directly the string order parame- 
ter in the Haldane phase with second order correlation 
functions, topological order can in principle be extracted 
from the higher moments of the output light quadrature. 
This measure, provided that the shot noise of the incom- 
ing light is small, employing for instance squeezed light, 
adds no experimental effort and will be the subject of 
future investigations. 
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